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~TITLE MTHSASIN § Aes Ww Point Arcsine routine | 
IDENT /1-008/ : File: MTHASIN.MAR Edit: JCw1008 


RADA ARRRRARRRASAARRAAAAASAAS ALAA RA ASSL EL ERE SARA R SERRE ERE S ER ER ESET SEES ESS 


COPYRIGHT (c) 1978, 1980, 1982, 1984 BY 
DIGITAL EQUIPMENT CORPORATION, MAYNARD, MASSACHUSETTS. 
ALL RIGHTS RESERVED. 


THIS SOFTWARE ia Sgt UNDER A LICENSE AND MAY BE USED AND gd 4 
ONLY IN ACCORDANCE WITH THE TERMS oF SUCH air tt air > ae pa, 3 


‘® ® 
:* * 
*® ® 
-* ® 
® ® 
:@ ® 
® ® 
*® ® 
* + 
ie RESON. NO TITLE TO AND OWNERSHIP OF THE SOFTWARE TS HEREBY = 
it TRANSFERRED. . 
** ® 
* THE INFORMATION IN THIS SOFTWARE IS SUBJECT TO CHANGE WITHOUT NOTICE  * 
** ® 
:® ® 
*® te 
** te 
:* ® 
*® ® 
*® & 
:® ® 


AND SHOULD NOT BE CONSTRUED AS A COMMITMENT BY DIGITAL EQUIPMENT 
CORPORATION. 


DIGITAL ASSUMES NO RESPONSIBILITY FOR THE USE OR RELIABILITY OF ITS 
SOFTWARE ON EQUIPMENT WHICH IS NOT SUPPLIED BY DIGITAL. 


LAR RRBRASRARARABAAAARARAAARALEARLALALESESESAR EARS RRR RRR EERE RAR ARR RRR RRR RR RAS SD SS 


: ABSTRACT: 


: MTHSASIN is a function which returns the floating point _arcsine in 
; radians of its single precision floating point argument. The call is 
; standard call-by-reference. 


; MTHSASIND is a function which returns the floating point arcsine in 
: ; defrees of its single precision floating point argument. The call is 
standard call-by-reference. 


; VERSION: 01 
: HISTORY: 
: ; AUTHOR: 
Peter Yuo, 29-Jun-77: Version 01 


; MODIFIED BY: 
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o v04-00 P 
; Detailed Current Edit History “SEP-19 are 3) 


; Flo r 
HISTO MTHRTL.SRCIMTHASIN.MAR; 1 
-SBTTL HISTORY ; Detailed Current Edit History 


] 

| 
ALGORITHHIC DIFFERENCES FROM FP=11/C ROUTINE: none | 
Edit History for Version 01 of MTHSASIN 


0-2 MTHSSERROR chenged to MTHS$SIGNAL. 
MTH$_... changed to MTH__.... 
Changed error handling mechanism. Put error result in RO before 
call ng MTHSSSIGNAL in order to allow user modify error Ad 
- Update version number and cone} ht notice. ; S 16-NOV-7 
Change MTH__INVARG to MTHSK JBS 9 -DEC-78 
Add "_"' to the PSECT directive. S 21-DEC-78 
Declare externals. SB, 17-May-19 
Use MTHSSQRT_RS5. SBL_ 27-Sept-1979 
Change JSB entry to MTHSASIN_R4. RBG 28-Sept-1979 
Added degree entry points. RNH 29-MAR-1981 
With the aid of the RTL error plotting routine is 
that MTHSASIN is more accurate for 'Xi>=0.65286 
evaluated by using (1-X)#(14X) instead of 1-x*2 
MTHSACOSD was found to be more accurate for :Xi 
1-X*2 was evaluated by (1-X)#(14x). 


8 ee tt 
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; Declarative Part of Modul MTHRTL.SRCJMTHASIN.MAR; 1 
-SBTTL DECLARATIONS ; Declarative Part of Module 


t INCLUDE FILES: 


EXTERNAL SYMBOLS: 


GBL 
sEXTRN MTHSSORT _R3 
TRN MTHSATAN- Re 
TRN MTHSATAND_R 
TRN ATHSK. SNVARGHAT 
TRN MTHSSSIGNAL 


EQUATED SYMBOLS: 


SF_1.0 = *F1.0 i 1.0 
LF_PI_OV_2 = *0007733a16 + 0) 
LF.90 = *x00004384 

value = 4 


: value.rf.r 


-PSECT _MTHSCODE PiC,SHR, LONG, EXE ,NOWRT 
7 program section for math routines 


OWN STORAGE: none 


CONSTANTS: 


tH HH HK HQOOODOOOCOCOCOOOOOOOOOOOWWMODGDOCDODCDODCDCO—“y *- 2 
NOUS WN (OOD NAU EWN SO OD NA UE WN 0 OONOU SEWN O0 
. 


00004080 
OF DB40C9 
00004384 
00000004 
: 
0 


tt 4 ss as = 2 - ss 


MACROS none 
PSECT DECLARATIONS: 


Ll 9 
MTHSASIN nt Arcsine routine 16-SEP-1984 4: AX/VMS Macro V04-00 Pa 4 
1-008 ndard Single Precision FL 6=SEP=1984 atk i 3} ERTHRTL. SRCIMTHASIN. MAR: 1 ” (4) 


-SBTTL MTHSASIN = Standard Single Precision Floating ASIN 


=z: 
—_ 
z=" 
Ana- 
ro 


++ 

FUNCTIONAL DESCRIPTION: 

; ASIN = single precision floating point function 
ASIN(X) is computed as: 


If X = Q, then ASIN(X) = Q. 

If X = 1, then ASIN(X) = PI/2, 

If X = oA. then ASIN(X) = =PI/2. 

If 0 < 1X! <1, then ASIN(X) = ATAN(X/SQRT(1-X##2)). 
If 1 < iXi, error. 


CALLING SEQUENCE: 
ASIN.wf.v = MTHSASIN(x.rf.r) 
; INPUT PARAMETERS: 


LONG = 4 3; define longword multiplier 
x = 1 * LONG : Contents of x is the argument 


IMPLICIT INPUTS: none 
OUTPUT PARAMETERS: 
VALUE: floating arcsine of the argument 
IMPLICIT OUTPUTS: none 
COMPLETION CODES: none 
SIDE EFFECTS: 
Signals: MTH$_INVARG if iX! > 1 with reserved operand in Re (copied to 
the signal mechanism vector CHFSL mi H_RO/R1 by LIBSSIGNAL 


Associated message is: "INVALID ARGUMENT’. Result is Foonein ed 
operand -0.0 unless a user supplied (or any) error handler changes CHFSL_MCH_RO/R1 


00000004 
00000004 


NOTE: This procedure disables floating point underflow, enables integer 
overflow. 


Be Se Ge Se Ge Ge Ge Ge Ge Ge Se Ge Ge Ge Ge Gee Ge Se Ge Ge 


OOOCOCSOOCSOOCOSOOCOOOSOSOSOSOSOOSOSOOSCOSOOOOOSOOOSOOSOOOOSOOOSOOOOOSOoOS = 


SOOoQOOoooooooo 


-ENTRY MTHSASIN, “M<IV, R2, R3, R4> 

3 standard call-by-reference entry 

; disable DV (and FU), enable IV 
MTHSFLAG_JACKET : flag that this is a jacket procedure in 


MOVAB G*MTHSSJACKET_HND, (FP) 


401C 


a a a a a ad a UD UO 
SOO NAN NN OS OO NAME WIN  O OOO NA UNE WN  O OD NAME WO VONOAUE WOO | $4 


i 
t 
| 
3 
3 
3 
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4 
4 
4 
4 
4 
4 
4 
4 
4 
5 
5 
5 
5 
5 
5 
5 
5 
5 
5 
6 
6 
6 
6 
6 
6 
6 
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set handler address to jacket 
handler 
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SOOCOCCOCSOOSOOOOOOCOOOOOOOOOSOOOOOSOOSOOSOOOOOOOOOOOOCOOSOOOOOOOSOO “Oo 


Sooooooofo 


tAresine routine.” ru 'OrSECT}SRS 94:96:9) PATENTS Becta Nee fOans1 P89" dy 


wis" lassie.” 


Le 


2 
| i j se of an error in roytine 
i : and res igna convert signal to user PC 
50 04 Bt $0 9 17% MOVF  avalue(AP), RO RO = }X! = avalue(AP) 
ert 0 D 175 MTHSASIN_RG ; call special ASIN routine 
04 sta 176 z return with result in RO 


sc hemeeaaeiceteieilicsteser: ; 
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! 4 -SBTTL MTHSASIN_R4 = Special ASIN routine 
! } 1 ; Special ASIN = used by the standard routine, direct JSB call. 
1 1 Q ; CALLING SEQUENCE: 
1 184 ; save anything needed in RO:R4 
1 185 ; MOVF ane ; input in RO 
010 1 $: JSB MTHSASIN_RG 
4 1 3 RSB 3 return with result in RO 
O10 188 
4 190 MTHSASIN_R4:: ; special ASIN routine 
1 191 MTHSASIN_RS:: ; Release 1 name 
54 50 30 01 13 Mov RO, R4 3 save X in R4 
46 13 «0013)—=—19 BEQL RETURN : Return ASIN(O) = 0 
01 194 
01 195 ; 
01 138 3; 0< ix! 
0015 197; 
0015 198 
50 8000 8F AA 0015 199 BICW #*x8000, RO 3; RO = ix! 
08 50 51 OO1A 200 CMPF =soaRO, S*#SF_1.0 + compare !X! with 1.0 
2c 18 001D 01 BGEQ GEQ_T0_1.0 3; branch if iXi >= 1.0 
OO1F O 
001F 03 ; 
OOF 04; 0< ixi ¢ 1.0 
SE 
22304027 8F 50 51 OO1F 07 CMPF RO gnygesesee? 3 compare ‘Xi with 0.6528662 
0C 19 0026 208 BLSS_ L$§_0.6 > branch if !X! < 0.6528662 
51 08 50 43 0028 09 SUBF3 RO, S*#SF_1.0, R1 ; RI=1-xX 
50 08 40 002c 210 ADDF2 S*#SF_1.07 RO : RO=1 +X 
50. 0 51 44 002F 11 MULF R1, RO ; RO = 1.0 = xX*2 
07 +11 «0032 1 BRB FLOW : join the normal flow 
bnee 135 LSS_0.65: 
50 50 44 0034 14 MULF RO, RO 3; RO = Xee2 
50 08 50 43 Boer 15 SUBF RO, S*#SF_1.0, RO 3 RO = 1.0 = X*#2 
038 16 FLOW: 
90009000" Er ig 038 17 JSB MTHSSQRT_R3 ; RO = SQRT(1-X*#2) 
54 0 4&4 041 18 DIVF3 P e 3: RO = X/SQRT(1-X*##2) 
00000000‘ EF 17 Onee i? JMP MTHSATAN_R4 3 RO = ATAN(X/SQRT(1-X**2)) 
0048 1 
pee § 3 
48 s 7 wé€ 78) 
pe 4; 
48 5 
pee § GEQ_T0_1.0: 
OF 14 re BGTR ERROR ; branch to ERROR if iXi > 1.0 
tb 599 
4D 0 3 
4D 13 ix! = 1.0 
tb 538 | 
50 OF DB40C9 ef 29 4D 4 MOVF #LF_PI_OvV_2, RO ; RO = PI/2 | 
4 54 5 TSTF = RG ; test the sign of xX | 
| 
} 
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1-008 rHsne ASIN "Re - Special ASIN rout ine emer rae £36; } ERTMRTL. SR SATHASIN.MAR; 1 ’ by 
1 BGEQ RETURN ; branch if xX > 0 
50 93 38 38 MNEGF RB RO ; RO = -PI/2 
3 RETURN: RSB 3; return with result in RO 
5 @ 
C r 3 1< Kt, error 
gE 548 
ge DD QO05C rf ERROR: PUSHL (SP) ; return PC from JSB routine 
7E 00' 8F 9A b25 MOVIBL #MTHSK —ess -(SP) ; 
50 01 OF 78 006 re ASHL #15, #T, RO ; RO = result = reserved operand -0.0 
B38 4 : goes to signal mechanism vector 
066 48 ; CCHFSL_MCH_RO/R1) so error handler 
066 49 : can modify the result. 
QOOO0000'GF 02 FB 0066 50 CALLS #2, G*MTHSS$SIGNAL : ON eb error and use real user's PC 
06D 51 : ependent of CALL vs J 
05 Bpee 26 RSB : rae - RO restored from CHFS$L_MCH_RO/R1 
006E 54 


| 
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-SBTTL MTHSASIND = Standard Single Precision Floating ASIND 


+4 
FUNCTIONAL DESCRIPTION: 


ASIND = single precision floating point function 


; ASIND(X) is computed as: 


If X = 0, then ASIND(X) = F 

If X = 1, then ASIND(X) = 90 

If X = ={. then ASIND(X) = -90. 

If 0 < 1X! < 1, then ASIND(X) = ATAND(X/SQRT(1-x##2)), 
If 1 < iki, error. 


; CALLING SEQUENCE: 


asind.wf.v = MTHSASIND(x.rf.r) 
INPUT PARAMETERS: 


\ 
LONG = 4 ; define Longuerd multiplier 
x = 1 * LONG ; Contents of x is the argument 


IMPLICIT INPUTS: none 
OUTPUT PARAMETERS: 
VALUE: floating arcsine of the argument 
IMPLICIT OUTPUTS: none 
COMPLETION CODES: none 
SIDE EFFECTS: 
Signals: MTHS_INVARG if iX! > 1 with reserved operand in RO (copied to 
the signal mechanism vector CHFS$L_MCH_RO/R1 by LIBSSIGNAL). 
Associated message is: ‘INVALID ARGUMENT'’. Result is reserved 
operand -0.0 unless a user supplied (or any) error handler changes CHF$L_MCH_RO/R1 


NOTE: This procedure disables floating point underflow, enables integer 
overflow. 


~-ENTRY MTHSASIND, “M<IV, R2, R3, R4> 

; standard call-by-reference entry 

; disable DV (and FU), enable IV 
MTHSFLAG_JACKET ; flag that this is a jacket procedure in 


MOVAB G*MTHSSJACKET_HND, (FP) ; 
set handler address to jacket 
handler 


MTHSASIN 
1-008 


ingle 


MOVF 
BSBB 
RET 


routine 
gle ine ston 4, 


@value(AP) 
MTHSASIND_ g,"° 


P=1986 $1 a 3 


AX/VMS Macro v04-00 Pa 
MTHRTL.SRCJMTHASIN.MAR; 1 i 


case of an error in routine 

: Tf an error, convert signal to user PC 
and resigna 

RO 2 ixis« @value (AP) 

call special ASIND routine 

return with result in RO 
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43 16 -SBTTL MTHSASIND_R4 = Special ASIND routine 
43 3 ; Special ASIND = used by the standard routine, direct JSB call. 
fa: 0 ; CALLING SEQUENCE: 
7E 1 3 save anything needed in RO:R4 
O7E § : OVF nt ; input in RO 
O7E ; JSB MTHSASIND_R4 
O7E 6; KSB 3 return with result in RO 
ee 
O7E MTHSASIND_R4:: 3; special ASIND routine 
54 50 $0 OO7E $28 MOVE RO, RG S seve Bin 
? Se | 008 ; BEQL D_RETURN : Return ASIND(O) = 0 
008 13 
008 3 70 < Ix! 
008 33 3 
008 ¢ 
50 8000 8F AA 008 5 BICW #*x8000, RO ; RO = ix! 
08 50 51 0088 6 CMPF =soaRO, « S*#S5F_1.0 > compare {Xi with 1.0 
2c 418 0088 7 BGEQ D_GEQ_T0_T.0 : branch if !X! >= 1.0 
008D 8 
008D 9; 
008D 40 ;0< ix! < 1.0 
008D 41 ; 
008D 4g 
B35C403A 8F 50 51 008D 4 CMPF RO, #*XB35C403A 3; compare {X! with 0.72929936 
0c 19 0094 44 BLSS L . 3; branch if ix! < 0.72929936 
51 08 50 43 0096 45 SUBF3) = RO. S*#SF_1.0, R1 sR1 es 1 - X 
50 08 40 009A 346 ADDF2 S*#SF_1.07 RO > RO=1+X 
50 51 44 009D 47 MULF R1, RO ; RO = 1.0 = x*2 
07 «11 0A) 48 BRB D_FLOW : join the normal flow 
OA 49 D_LSS_0.729 
50 50 44 OOA 50 ULF RO, RO 7; RO = Xe*2 
50 08 50 43 OOA 51 SUBF RO, S*#SF_1.0, RO 7 RO = 1.0 = Xee2 
00A9 3g D_FLOW 
OOO00000'EF 16 BoA? 5 JSB MTHSSQRT_R3 ; RO = SQRT(1-xX*#2) 
54 50 47 OOAF 54 DIVF3 ~ RG, 3 RO = X/SQRT(1-X*##2) 
00900000" EF 17 4 2? MP MTHSATAND_R4 ; RO = ATAND(X/SQRT(1-X**#2)) 
aid 2$ 
0B9 58 ; 
94 59; 1 =< Ix! 
9 60 ; 
089 61 
94 6 D_GEQ_T0.1.0: 
Al 14 9 6 BGTR ERROR ; branch to ERROR if ix! > 1.0 
088 64 
088 65 
8B 06 ; 
67 ; ix! = 1.0 
8 $s ‘ 
50 00004384 ef 39 % MOVF #LF_90, RO ;RO=90 | 
~ -% BC 71 TSTF RG 3; test the sign of X 
95 18 C4 72 BGEQ RETURN ; branch if xX > 0 
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nose 


atin Poin Arcsine rout 
InDoRG Special ASIND 


73 MNEGF 
73 D_RETURN: 


5 RSB 
7 


15 ~ END 


ine” 
routine 


RO, RO 


"S-SEp 
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; RO = -90 
3; return with result in RO 


Page 


MTHSASIN 
Symbol table 


LSs_0O. ‘ 
MTHSSJACKET_HND eeererer 
MTH$$S1IGNAL 


geereere 


MTHSASIN 9990000 
09 O0006E 
00000 

MTHSASIN_RG 90000016 

MTHSASIN-RS 0000001 

MTHSATAND RG etereres 

MTHSATAN R4 eeerereee 

MTHSK_INJARGMAT eeeeeces 

MTHSSORT_R3 eereeeese 

RETURN 00000058 

1.0 = 00004080 

vaCuE = 00000004 

PSECT rame 

. ABS . 

_MTHSCODE 

Phase 

Initialization 


Command processing 
Pass 

Symbol iable sort 
Pass 2 

Symbol table output 
Psect synopsis output 


Cross-reference output 
Assembler run totals 


The workin 
5683 ° 
There were 
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. 
x OO 
G 01 
RG 4 
RG 1 
G 01 
RG 01 
x 00 
x 00 
x 00 
x 00 
i 01 
er ee een eece see eet 
! Psect synopsis ! 
gem e now en enn ane cne + 
Allocation PSECT No. Attributes 
00000900 g-} 00 ¢ 0.) NOPIC USR CON ABS ~ LCL NOSHR NOEXE NORD 
OQOOOOOCA (¢ 202.) O01 ¢ 1.) PIC USR CON REL LCL SHR’ EXE RD 
tee eee ret ee eee ewe cee eee $e 
H Performance indicators ! 
Page faults CPU Time Elapsed Time 
00:00:00.09 0:00:00.74 
133 0:00:0 99 0: 8:8 89 
9 :00:00.97 3: 8: .61 
7:00:00. :00:00.07 
7 :00: $ 0:00:02.6 
:00:00.04 0:00:00.1 
5:00; 0:00 §:00:00:00 
331 00; :02. 3 00; 0:11.19 


st Limit was 750 pages. 
bytes (12 pages) of virtual memory were used to buffer the intermediate gepe. 
10 pages of symbol table space allocated to hold 26 non-local and 0 local symbols. 
439 source Lines were read in Pass 1, producing 14 object records in Pass 2. 
1 page of virtual memory was used to define 1 macro. 
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ASIN.MAR;1 


NOWRT NOVEC BYTE 
NOWRT NOVEC LONG 


~00 Page 1 
eS (5). 


ns 
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Br wes in 
VAX-11 Macro Run Statistics MTHRTL.SRCJMTHASIN.MAR; 1 (7) 


oeeeeceeceeeeeeaeeseesseseeese + 


! Macro Library statistics ! 


Goeeeoeecoeeeoeoceeecoeeseeceocacees} 


Macro Library name Macros defined 
~S255SDUA28: (SYSLIBISTARLET.MLB; 2 See Ke 
0 GETS were required to define 0 macros. 

There were no errors, warnings or information messages. 


MACRO/ENABLE=SUPPRESS1ON/DI SABLE=(GLOBAL , TRACEBACK) /LIS=LIS$:MTHASIN/OBJ=OBJ$:MTHASIN MSRC$:MTHJACKET/UPDATE=(ENHS:MTHJACKET) *MSRCS: 
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